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Abstract 

We introduce a formalism for the calculation of the time of arrival t at a detector of 
particles traveling through interacting environments. We develop a general formula- 
tion that employs quantum canonical transformations from the free to the interacting 
cases to compute t. We interpret our results for the time of arrival operator in terms 
of a Positive Operator Valued Measure. We then compute the probability distribu- 
tion in the times of arrival at a detector for those particles that - after their initial 
preparation - have undergone quantum tunneling or reflection due to the presence 
of potential barriers. We make a extensive analysis of several cases, the main results 
being the presence of the expected retardation or advancement for transmission, and 
of non-foreseen two bump structures for some cases of reflection. 

03 
US 

cr. 

> 



On 

S3 

Oh 



IMAFF-PCC/99-17 
March 1999 



iWork supported in part by DGESIC PB97-1256. 

2 Instituto de Ffsica, UFRJ. 21945-970 Rio de Janeiro- Brasil. 

3 Departamento de Fisica, Universidad de Alcala, Alcala de Henares, Spain, 
e-mail addresses: julve@fresno.csic.es, leon@ifF.csic.es, pitanga@if.ufrj.br, 
fcrnando.urries@alcala.es 



1 Introduction 



In this paper we work out a theoretical framework to compute the time in which a particle 
that moves in an interacting environment arrives at a given point. We apply our results 
to answer the long pending question: How long does it take for a particle to tunnel across 
(or through) a quantum barrier? In the construction of this framework we will have to 
deal with problems of very different kind that we introduce now: 

First, there is the nature of time in quantum mechanics. It appears as the external 
evolution parameter in the Schrodinger and Heisenberg equations, common to both, sys- 
tems and observers alike. However, time arises in many instances (transitions, decays, 
arrivals, etc.) as a property of the physical systems. The attempts to promote time to the 
category of observable run early into the obstruction detected by Pauli [Q: A self adjoint 
time operator implies an unbounded energy spectrum. This was soon related to the uncer- 
tainty relation for time and energy, whose status and physical meaning has been subject 
of controversies 0, [| |], |5|, and is still subject of elucidation today (see for instance || 
and f7j). The question remains unsettled for closed quantum systems, specially in the case 
of quantum gravity, whose formulation is pervaded by the so called problem of time ||. 

Second, the definition of the time-of-arrival (toa), which is probably the simplest can- 
didate time to become a property of the (arriving) physical system, rather than a mere 
external parameter. Due to its conceptual simplicity, it has been used in many cases to 
illustrate different problems related to the role of time in quantum theory. Allcock ana- 
lyzed || extensively the difficulties met by the toa concluding they were insurmountable. 
The present situation is ambiguous. On the one side, there are sound theoretical analy- 
sis (TO, ry| of the toa showing that it can not be precisely defined and measured in quantum 



mechanics. This contradicts the possibility of devising high efficiency absorbers [[HJ, that 
could be used as almost ideal detectors for toa [I3|. On the other hand, there are ex- 



plicit constructions of a self-adjoint toa operator for the non-relativistic free particle in 
one space dimension []H , the relativistic free particle in 3-D both avoiding the Pauli 



problem. There is also an alternative formulation |16| as a Positive Operator Valued Mea- 



sure (POV). Finally, the toa has been measured in high precision experiments [17, 
on the arrival of two entangled photons produced by parametric down-conversion, one of 
which has undergone tunneling through a photonic band gap (PBG). The experimental 
results that show superluminal tunneling, neatly identify the Hartman effect |19j and the 
Wigner time delay [^(J (or phase time) as the physically relevant mechanisms for the tun- 
neling time and toa respectively. Whether these results apply only to photons and are 
due to the specific properties of the PBG used, or can be extended to other particles and 
barriers, can not be decided in the lack of a satisfactory theory of the toa at a space point 
through interacting environments. 

The third question is thus the tunneling time, for which there are three main proposals. 
Wigner introduced the phase time in his analysis |2(J of the relationship between retar- 
dation, interaction range, and scattering phase shifts. Buttiker and Landauer introduced 
the traversal time in their study of tunneling through a time- dependent barrier. Soon 
after, Buttiker used the Larmor precession as a clock | 22fl , identifying the dwell [p3f| , traver- 
sal, and reflection times as three characteristic times describing the interaction of particles 
with a barrier. Recent reviews that include these and other approaches, discussing toa and 



1 



tunneling times from a modern, unified perspective, can be found in [23] and |p5 |. The 
light shed on these questions by the two photon experiments is revised in [^] and [ffifl. 

The plan of the work is as follows: In Section 2 we construct the generic toa formalism, 
generally valid, but still awaiting the necessary connections to each physical situation 
of interest. The starting point is the case of the free particle. A suitable canonical 
transformation, similar to that used in scattering theory, gives the toa in interacting 
environments, (also at points where V(x) 7^ 0). In Section 3 we specialize to the case 
of potential barriers, showing the existence of two characteristic times of arrival at each 
point, one for incoming and other for outgoing particles, in the same way as there are two 
different kinds of scattering states associated to the free "in" and "out" asymptotic states. 
In Section 4 we make some calculations and give estimates of interest for the case of the 
square barrier. Here we give the toa for both, transmitted and reflected particles, and 
discuss some properties of the associated probability distributions. Finally, we summarize 
our results in Section 5. 



2 Time of arrival formalism 

To measure the time of arrival of a free particle at a point x one would: a) place a detector 
at x, b) prepare the initial state of the particle at t — 0, and then, c) record with a 
clock the time t when the detector clicks. The value of t gives the toa of the state 
at x. Repeating this procedure with identically prepared initial states, one would get the 
probability distribution in times of arrival at x. Of course, the results would depend on 
the initial state chosen, which stores all the information regarding the initial distribution 
in positions and momenta of the particle. 

To determine the effect on these times of climbing (or tunneling through) a potential 
barrier, one would simply put the barrier in between the detector and the initial state, and 
then record the new times of arrival. With an initial state identical to that prepared for the 
free case, any difference in the probability distributions should be an effect of the barrier. 
Several questions can be investigated by changing the properties of the barrier: its height 
or width if it is rectangular, even its very form. This has been explicitly done in the two 
photon experiments at Berkeley, by putting alternatively a mirror and an ordinary glass 
in the path of one of the photons. It would also be of interest to analyze the dependence 
on x when V(x) > 0, i.e. with the detector within the range of the interaction. 

In classical mechanics particles move along the trajectories H(q,p) =const. as t in- 
creases. This allows to work out t x , the time of arrival at the point q(t) = x, by identifying 
the point (q,p) of phase space where the particle is at (say) t = 0, and then following the 
trajectory that passes by it, up to the arrival at x. The mathematical translation of this 
procedure is given by the equation of time: 

rn f x dq' 



t x (q,p) = sign(p)^j — J 



1 ^H(q,p)-V(q>) 

that is discussed at length in many textbooks. 

For free particles Eq. (|Tj) gives t®(q,p) = ra[x — q)/p that, in spite of its seemingly 



simple and harmless appearance, presents some problems in its quantization [EL 14, [L5|. 
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First of all, it requires symmetrization: 



i°MP) = m(- - ±{q, -}+) = -e-^MqMe^, (2) 
p 2 p y p y p 

As is well known, the eigenstates of this operator \txsO) in the momentum and energy 
representations can be given as (h = 1) 

I p p 2 

(pltxsO) = 6(sp)\ exp(i — t — ipx) (3) 

V 2nm 2m 

where we use s = r for right-movers (p > 0), and s = I for left movers (p < 0.) The 
indexlabel stands for free case. Finally, the argument sp of the step function that appears 
in the momentum representation is +p for s = r, and — p for s = I. The degeneracy of the 
energy with respect to the sign of the moment is explicitly shown by means of the label s 
in the energy representation 

(Es'0\txs0) = 5 s , s -^= exp(iEt - %s\/2mEx) (4) 
V 2n 

where the s in the exponent of the rhs stands for +1 for s = r and —1 for s = I. These 
eigenstates are not orthogonal, which in the past gave rise to serious doubts about their 
physical meaning. The origin of this problem can be traced back to the fact that (fj) is 
not self-adjoint, that is {ip\f^}) ^ (t®ip\ip). This was proved by Pauli long time ago 
and is due to the lower bound on the energy spectrum. The problem emerges as soon as 
one attempts integration by parts in the energy representation. Ref. |25] contains a recent 
illuminating review of these and other related questions. 

A close inspection of (^]) reveals that it still needs of some prescription to deal with 
p = 0, in order to produce a self-adjoint operator with orthogonal eigenfunctions. When 
p — > + the t° x eigenstates span the subspace {\Er)}, that comes abruptly to an end at 
p = 0, and lacks the continuation to negative values of p necessary to get orthogonality 
(in other words, a delta function coming from the integration over p would require p G 



(— oo, +oo), but actually p G (0 + , +oo) alone). The way out [|Tj] is to stretch analytically 
the interval p = (0, e) to p = (— oo, e) in such a way that the new values of p cover all the 
real line. The same problem, with the same solution, occurs when p — ► 0~. Observe that 
the attempt to "glue" {\E r)} to {\E I)} in a neigbourhood of p = is misleading, because 
it would connect with continuity states belonging to macroscopically different experimental 
setups (those moving from the far left to the far right with small but positive momentum, 
with those moving from the far right to the far left with small but negative momentum). 
In fact, there is no discontinuity in the eigenstates \p) at p = while, at the same time, 
{\E r)} can not be connected continuously to {\E I)}. This is a consequence of the different 
physical meaning of these two classes of states, given mathematically by different kinds of 
boundary conditions. The correspondence principle maps classical trajectories to energy 
eigenstates, but does not restrict the momentum eigenstates in the same way. 

The measurement problem associated to the toa operator can be solved more simply by 
the use of a Positive Operator Valued Measure (POV) , that only requires the hermiticity 
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of i° x (i.e. t° = (t°)* T ). Here, instead of a Projector Valued spectral decomposition of the 
identity operator, one has the POV 

r 2 

Px(ti,t 2 )=Y dt\txsO)(txsO\ (5) 

where P(l,2) 2 7^ P(l,2) because \txsO)(txsO\ is not a projector, as the states are not 
orthogonal, but where the limit as t — > 00 of P(—t, +t) is the identity. The attained time 
operator is no longer sharp, but is well suited for measurement. This solution has been 



implemented in |L6| , and extensively analyzed in refs. |2q , |29fl and in the review 

In general, Eq.([l|) can not be explicitly solved in the presence of interaction. This 
is not a problem classically, but makes hopeless any attempt to straight quantization. 
Even for the simple rectangular barrier that looks similar to a piece-wise assembly of 
free cases, there is no successful way to employ ([I]) directly. Here, the reason is not the 
algebraic complexity of that expression, but the fact that it only applies for values of x 
that are classically within the reach of (q,p), while in quantum mechanics all values of 
x are possible. Thus, neither quantum tunneling nor quantum reflection have a classical 
analog in (ffl), (that gives complex numbers for these cases). Lacking it, this equation is 
a troublesome starting point for quantization. In spite of these problems, it arises as a 
critical time for the behaviour of particles in time dependent barriers, in the form of the 
so-called traversal time tt'- the time a particle takes to traverse the barrier |21[ (which 
is the absolute value of ([!])). As was indicated in the introduction, this time has been 
compared, giving interesting information about tunneling and related phenomena, 



with the "dwell time" [£3] r d , defined as the ratio of the number of particles inside the 



barrier to the incident flux, and the "phase time" [2(| r^, the time associated to the peak 
of the propagating wave packet by the stationary phase approximation. 

In this work we desist from the attempt of straight quantization of the classical expres- 
sion. Instead, we will construct the solution to the interacting case in terms of the well 
known results that apply to the free case. The aim is similar to what has been done in for- 
mal scattering theory, where the Lippmann-Schwinger equations give the complete Green 
function and the scattering amplitude in terms of the free Green function and the poten- 
tial. This procedure also serves to get the energy eigenstates (stationary scattering states) 
in terms of plane (or spherical) incoming or outgoing waves. Of particular interest for our 
problem is the existence of the Moller operators Q± that connect the free particle Hilbert 
spaces TCi n a.ndH out , to the Hilbert space Tt of the scattering and bound states. These 
operators are only isometric in the presence of bound states, because the correspondence 
between states in TC and free states can not be one to one. In this paper we will consider 
only well behaved potentials (V(q) > Vg, V(q) vanishes at the spatial infinity, and has no 
valleys to avoid resonances), for which the operators are unitary because there is one free 
state for each scattering state. In this case, the intertwining relations HQ± = Q±H can 
be put in the form H = Q±H Q±. This is an useless relation as it stands, but it prompts 
the following question: Is there a similar relation for obtaining the time of arrivaltoa in 
the presence of interactions in terms of i°? Below, we construct a positive answer, giving 
explicitly the generic solution for the positive potentials we are considering. In the next 
section we will apply our results to different cases of interest. There, we will make the 
appropriate distinctions between incoming (+) and outgoing (— ) cases, that are necessary 
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to predict the probability distributions, but that do not affect the toa, or its eigenstates, 
when there is a unique free particle Hilbert space, i.e. when TL in = H ou t as in our case. 

We now come back to phase space and recall that the classical canonical transforma- 
tions q = q(q,p), p = p(q,p) can be defined implicitly by the use of auxiliary functions 
F, G, F, G in the following way: 

F(q,p) = F(q,p) (6) 
G(q,p) = G(q,p). (7) 

It is easy to work out the following relation among Poisson brackets: 

{F > G}g p{q,p}q p = {F,G} q p (8) 

In these conditions, the transformation is canonical (i.e. {q, p} q p — 1) if and only if 

{F , G}q p = {F,G} g p . (9) 

This relation has the additional property of fixing one of the four functions F, G, F, G, 
once the other three are given. We can choose F and G as the free particle Hamiltonian 
and time of arrival respectively. Then, if F is the complete Hamiltonian H, G will be the 
corresponding toa t x , and the above equation gives:{t x , H} = 1, whose solution along the 
classical trajectories is (p. 

Canonical transformations were introduced by Dirac in quantum mechanics by the use 
of unitary transformations U (UU^ = U^U — 1 ). If the operators q,p are canonically 
transformed from q, p , then there is a unitary transformation U such that 

q = UqU ] 

p = UpUl (10) 
Then we can define implicitly quantum canonical transformations, like the classical ones, 

byipj 

F{q~p) = UF{q,p)U ] = F{q,p) 

G(l~p) = UG(q,p)U^ =G(q,p) (11) 

where the last equality in each row is the definition of the barred operators in terms of 
the unbarred ones, while the first equality comes from the straight application of fllOf) to 
the l.h.s. Being U a unitary transformation, the spectra of the canonically transformed 
operators have to coincide, that is: 

a(q) = a{q) = K, a(p) = a(p) = K 

a(F) = a(F), a{d) = a(G) (12) 

where the second row stands because F and F, G and G are also unitarily related operators. 
We will now show how to use the above relations to build the operator G once F, G and 
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F are given. The only restriction to our solution is that the operators F and F are self- 
adjoint, so that their eigenstates are orthogonal and form a complete basis. The eigenstates 
of these operators corresponding to the same eigenvalue Xf are: 

F\f)=X f \f), F\f) = X f \f) (13) 

They form orthogonal and complete bases, i.e. they satisfy 

(fs\f's')= 5 ss ,5(X f -\> f ),T,sUx)d\f\fs)(fs\ =1 (14) 
(fs\f's')= 6 ss ,6(X f -X f f ), EsL { x)dX f \fs)(fs\ =1 (15) 

where we allow for some degeneracy (that has to be the same for both F and F) labeled 
by s. We have also assumed that A is continuous, while s is a discrete index. These 
assumptions could be changed straightforwardly if it were necessary. Now, an operator U 
satisfying the first row of Eq . flTTD can be given simply as: 



tf = E /, d\ f \fs)(f 8 \ (16) 

It is straightforward to verify that it is unitary. We can now proceed to the sought for 
result: the definition of G in terms of G using U, that is G = WGU. The full fledged 
expression is 

G(q,p) = £ / dX f dX r \fs) (fs\G(q,p)\fs>) (f's'\ (17) 

that constitutes our main result in the quantum canonical formalism. We will now apply it 
to the case where F is the free Hamiltonian H , F the complete Hamiltonian H: F(q,p) = 
p 2 /2m + V(q), and G the time of arrival of the free particle Eq.([|). Then we have 

U (±) = J2 dE\EsO)(Es(±)\ (18) 

where the label in \EsO) indicates that it is an eigenstate of H , its absence corresponding 
to the eigenstates \Es) of the complete Hamiltonian. For notation simplicity, we defer to 
Section 4 (Eq. (|58l )) any reference to the incoming or outgoing (±) character of the time 
operator and states we are constructing, as it is unnecessary here. After substitution in 
Eq.([l7]) we get for the toa at x in the presence of the potential V 

POO 

t x (q,P) = Y dEdE'\Es) (EsO\i° x (q,p)\E's'0) (E's'\ (19) 

where t° x is the toa for the free problem, as given by @. The bracket in the above 
expression involves free quantities only and can be computed to give 



(EsO\t x=0 (q,p)\E's'0) = -(EsO\^fq^f\E's'0) 

= -^ss> (&{E - E')) (20) 
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where we have used (p\EsO) = y^<5(p — sy/2mE), and put x = for simplicity (we will 
restore the x dependence later on). Using the above expression in (|19|) , we get the operator 
form 

«-«.» = ~Ef l*> (^) W < 21 > 

that allows the explicit evaluation of the toa for all the cases in which the solutions 
^Esio) — (qIEs) to the Schrodinger equation are known. 

By solving the eigenvalue problem posed by (0), we would get the spectrum and 
eigenf unctions of the toa. Instead of using that procedure, we will do it in all generality 
(within our assumptions about the potential and Hilbert spaces) by using the canonical 
transformation U to get the answer from the spectrum and eigenf unctions of the free case. 
From the positive operator valued decomposition of Eq. @ we can write 

/+oo 
dtt\tsO)(tsO\ (22) 
-co 

where we still keep x = 0, and the spectrum of i° x is R. The states \ts0) are the non 
orthogonal eigenstates of this operator given in (£|). When used in (|20|), Eq. (p^) leads to 

/+oo 
dtt(EsO\ts"0)(ts"0\E's'0) 
-co 

1 r+oo 

= 7^8 SS > dtt exp(i(E - E')t) (23) 
2ix J-oo 

Putting this expression in Eq.fliTS|) we get the spectral decomposition for the toa at x = 
in the presence of the potential V(q): 

tx=o(q,p) = J2 dtt\ts)(ts\, (24) 

„ J —oo 



where the toa eigenstates are defined by 

1 r+oo 

These eigenstates are not orthogonal. By direct computation using ( |25| ) one gets 

1 i S sx i 



1 r+oo 

Its) = —= \ dE exp(iEt)\Es). (25) 
\/2tt Jo 



I r+oo 

(ts\t's') = —5 SS , / dE exp(-iE(t - t' - ie)) 
2tt Jo 



2tt f - t + ie 

/+co r+oo 
dt\ts)(ts\ = J2 dE\Es)(Es\ = I (26) 
-co „ Jo 



A set of relations defining the same kind of positive operator valued measure that appears 
in the free case [|16| . We can now restore the dependence of the different quantities on the 
detector position x. Restarting from Eq. (p0[) , but this time keeping the term mx/p in the 
brackets, we arrive to 

/ + CO 
dtt\txs)(txs\, (27) 
-co 
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with the position dependent eigenstates given by 

/-+oo 
_.. 



1 r+°° , 

\txs) = / dE exp (iEt-isxV 2mE)\Es). (2* 

\/2n Jo 



where the s in the exponent of the rhs stands for +1 when s = r, or —1 when s = I. We 
are now provided with the tools necessary for a physical interpretation. Given an arbitrary 
state ip at t — 0, its time of arrival at a position x has to be, according to (JZTD, 

/+oo 
dtt\(txsm 2 , (29) 

with the standard interpretation of 

£|<te#>| 2 (30) 

s 

as the probability density that the state arrives at the detector placed at x in a time 
t. 



3 Classical trajectories and energy eigenstates 

In three space dimensions, the free particle Hamiltonian Ho is independent of the posi- 
tion and of the direction of motion as a consequence of the homogeneity and isotropy of 
space. In the one dimensional case we are considering in this paper, Ho continues to be 
independent of the position, while the only remainder of isotropy is the degeneracy of the 
Hamiltonian in the sign of the momentum. For a given energy, the particle can move to 
the right or to the left, and occupy any position in space. With the passage of time (from 
— oo to +oo), the particle will successively pass by all the positions in configuration space. 
A steady flow of (say) one particle of energy E per unit time coming from the left - and 
necessarily going to the right at the same rate - will correspond to a properly normalized 
right mover stationary state. This applies classically in phase space, as well as quantum 
mechanically in the Hilbert space, and is also valid for left movers. 

In the presence of interaction, given by a potential energy V(q), space is no longer 
homogeneous, and the removal of degeneracies is more involved. This is the case for the 
potential barriers we are interested in. To simplify the analysis, avoiding bound states 
and resonances whose effect we study elsewhere, we will only consider here barriers with 
positive potential that vanishes sufficiently fast at the spatial infinity. In addition, these 
barriers will have at most a maximum of value V, and no minima. Two textbook barriers 
of this kind are the square barrier V(q) = V if < q < a, else V(q) = 0, and the well 
behaved barrier V(q) = V/ cosh 2 (q/d), whose energy spectrum and eigenstates are well 
known. They can be used to exemplify the material discussed in this section. 

The classical trajectories associated to these barriers are simple to compute, they are 
sketched in the phase space plane (q,p) in Fig.l, where the arrows on the trajectories 
indicate the direction of motion. A main feature clear in both, the upper and lower 
drawings, is the symmetry (q,p) «-> (q, —p), which comes from the degeneracy in sign(p). 
The symmetry under (q,p) <-> (—q,p) is due to the spatial symmetry of the potential 
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chosen, and is largely irrelevant for the present analysis. Another remarkable feature is 
given by the two symmetric curves crossing the origin. Consider first in the upper figure 
the curve that we call py of equation p = pv(q), where pv(q) =sign(g)y / 2m(V — V(q)). 
It divides the plane into two regions. All the trajectories incoming to the barrier from 
the left, and only them, are above p v . These trajectories share the property that as 
t — > — oo, (q,p) — > (—00, \/2mE), and split into transmitted trajectories in the case when 
E > V and reflected trajectories when E < V. Below pv, we have the physically equivalent 
set of the trajectories incoming to the barrier from the right and only them, verifying that 
as t — > —00, (q,p) — > (00, —\j2mE\ As above, there are also transmission or reflection in 
this region. The important thing is that the curve p = pv(q) serves to separate the phase 
space into two distinct (physically equivalent) regions in terms of the asymptotic initial 
condition at t — > —00. That is, above p v there are incoming right-movers and incoming 
left-movers below it. 

We could also make a partition of the phase space as shown in the lower part of Fig.l, 
using the curve p = —pv(q). Above it, we find all the outgoing right-movers; below it, all 
the outgoing left-movers. There is also transmission or reflection depending on whether 
E > V or not. The partition of phase space is here in terms of the final asymptotic 
condition for t — > +00: above — py there are right-movers ((q,p) — > (00, y2rnE)), and 
left-movers ((q,p) — > (—00, —\ / 2mE)) below it. Finally, we can count particle trajectories 
in two alternative ways, as incoming or as outgoing trajectories, according to the value of 
the asymptotic momentum at t — > — 00 or at t — > +00 respectively. For each way there 
will be right and left movers according to the asymptotic sign of p. The usefulness of a 
specific counting procedure will be dictated by the objective pursued. For instance, to 
investigate properties of the outgoing particles collected to the right (left) of the barrier, 
we have to select from the set of trajectories above (below) — p v , a selection in terms of 
Pv being useless to this end. 

In quantum mechanics one finds energy eigenstates in the Hilbert space instead of the 
classical trajectories in phase space. A particle will be described by a wave packet ip, 
a suitable superposition of energy eigenstates, that evolves according to the Schrodinger 
equation. If at some instant t — the state of the particle is \ip), the state at an arbitrary 
time can be written, using completeness of the energy eigenstates, as 

|^(0> =exp(-ifft)l^) =J2 I dEexp(-iEt)\Es(±)}(Es(±)\ilj) (31) 

s J 

where s — r,l gives the (asymptotic) sign of the momentum, necessary to remove the 
degeneracy of the energy eigenstates, that appears because E is independent of s. The 
signs (+) or (— ) indicate that the decomposition of the identity has been made in terms 
of the corresponding energy eigenstates: \Es(±)) = fl± \EsO). These are the quantum 
analogs of the classical trajectories. 

We now analyze in more detail the two mathematically equivalent ways for counting 
the energy eigenstates. As we said, the (+) states represent a steady flow of particles going 
into the barrier, corresponding to the partition of classical trajectories according to py. 
These states can be expanded in terms of the energy eigenstates \Es(+)). The label s — r,l 
would correspond to the sign of the momentum in the free case. Here, the states with 
s = r (s = I) are initial right-movers (left-movers). Observe that in the presence of the 
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potential the condition E > V no longer guarantees transmission, nor E < V reflection. 
The \E s (+)) states are widely used tools for studying the effects of quantum barriers. 
The usual procedure is to prepare an incoming right mover (q\Er 0) = e wq , and determine 
the coefficients of transmission and reflection with the help of the matching conditions. 
For each particle in the scattering state \Er(+)), these coefficients give its probability 
amplitude of finally being transmitted to the other side of the barrier, or alternatively, of 
being reflected, bouncing back to the left. Consider the case of a barrier contained in the 
range (0,a). The (+) states can be written in configuration space representation as: 

(q\Er(+)) = 6(-q)(e ipq + R(p)e~ ipq ) + 6(q)6(a - q) Xr (p, q) + % - a)T(p)e ipq (32) 
(q\El(+)) = d{-q)T{p)e-™ + 9{q)9{a - q) X i(p, q) + % - a)(e*» + fl(p)e^)(33) 

where p = \ / 2mE, the Xr,i solve the Schrodinger equation in the barrier region, and R 
and T are given by the matching conditions at q = and q = a. A (not displayed) overall 
normalization factor \jraj 2np will insure a steady incoming flux of one particle per unit 
time. 

The alternative way of counting energy eigenstates is in terms of the (— ) states. These 
are the energy eigenstates that emerge from the interaction region towards the spatial 
infinity. The \Es(—)) are states going out from the barrier towards the right (r), or 
towards the left (I). They correspond to the classical partition according to —pv- They 
describe a post selection: that of the states that will be recorded with certainty by the 
appropriate detectors placed to the right (r case), or to the left (I case), of the barrier. 
By normalizing like the (+) states, there will be one record, i.e. one outgoing particle per 
unit time in the steady flow corresponding to each of these states. We can give them with 
the same notation as before: 

(q\Er(-)) = d(-q)T{p)e ipq + 6(q)d{a - q) Xr (p, q) + % - a){e ipq + R{p)e- ipq ) (34) 
(q\El{-)) = 6{-q){e- ipq + R{p)e ipq ) + 6{q)9{a - q) X i{p, q) + % - a)T(p) e -*«(35) 

We said above that the (+) and (— ) states give mathematically equivalent decom- 
positions of the unity E in terms of the energy eigenstates and their degeneracies. In 




that was used in (pT|). As already noted, the physical meaning of these states is very 
different. It is well known the form in which the (+) states are used in three space 
dimensions to describe the stationary scattering states, and to get the scattering cross 
sections from its asymptotic form at \q\ — ► oo. There, (q\p(+)) ~ e m + f(p,Q)— and 

(da/dQ) = \f(p,Q)\ 2 , with p the asymptotic incoming momentum and Q the angle (p,q). 
In one dimension, the direction given by Q reduces to r or I, and the scattering amplitude 
f(p, O) to the appropriate transmission T and reflection R coefficients. Less used, but 
also well known is the association of the (— ) states to a post selection of the final states. 
For instance, \p(—)) is the state that comes out from the interaction region becoming 
asymptotically (at t — > oo) the free state \p). There are several different cases possible. 
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They can be described by means of projectors, a formulation that we will use later on. 
Define 

U(Es(+)) = \Es(+))(Es(+)\ (37) 

Il(Es(-)) = \Es(-))(Es(-)\ (38) 

the projector U(Es(+)) selects the states approaching the barrier from the left (s = r), 
or from the right (s = I), with energy E. In the same form, U(Es(—)) selects the states 
coming out from the right of the barrier (s = r) , or from the left (s = I), with energy E. 

Let V be the maximum value of the potential energy. Classically, there will be full 
transmission through the barrier for E > V, and complete reflection for E < V. As already 
said, this is not so in the quantum case in which partial reflection and transmission will 
take place for values of k at which they can not occur classically. In general, the wave 
packets will occupy regions of space that are classically forbidden, giving rise to pure 
quantum phenomena like tunneling, etc. It is the use of the above projectors that will give 
us the control of these situations in the Hilbert space, so that we can properly address the 
time of arrival at a space point in the presence of barriers, mirrors, etc. 



4 Gaussian wave packets on square barriers 

This section is devoted to the theoretical predictions provided by the previous formalism 
for the results of the experiments with potential barriers. We will analyze the case of 
transmission through a barrier in the first place, and then will turn our attention to the 
case of reflection, as both are aspects of the same problem. To fix notations, we will 
assume in what follows a square barrier of height V and width a, in the interval (0, a), 
so that (|32| , |33| , [34]) , and fl35|) are the energy eigenstates in the coordinate representation. 
We also assume that the detector is placed to the right of the barrier (that is, x > a), 
and that the initial state at t — is given by a Gaussian wave packet of width Aq = 5, 
centered at go < 0, with mean momentum p Q > 0. The wave packet in configuration space 
and its Fourier transform are: 

= ( _L )(i/4) e -^ e -(^-^o) 2 (39) 

2710 

(p\^) = (^)(l/4) e -«5 2 (p- P o) 2 -*WO (=-0(p)) 



For appropriate values of qo,Po and 5, such that p 5 » 1 and 5 « \qo\, almost all the 
packet is initially at the left of the origin and moving with positive momentum towards 
the barrier. We use this simplifying assumption in our qualitative arguments, and in the 
intuitive description of the process given below (indicated in the formulas by the use of 
~ instead of =), while working with the full expression ( |3TID wherever necessary in the 
numerical calculations. 

Our experiment differs from the usual ones in that it involves a post selection. Usually, 
one prepares a well defined initial state and observes its distribution into the possible final 
states. This gives the probabilities of these outcomes. Here, we select the final state of 
the particles coming out from the barrier to the detector. The procedure is to collect 
all the particles in this state, registering their times of arrival. From these records we 
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construct the conditional probability [31] in time we are interested in (i.e. in the toa's of 
all the particles gathered at the detector and only of them). Given the state that of 
a particle placed at the left of the barrier at t = as in (^), it is only the component 
= J dE U(Er(— )) \ip) that will come out to the right of the barrier (c.f. Eq. (PB|) ) . 
arriving at x > a in some time t. Our formalism gives for the probability distribution in 
these times of arrival 

P{-)(t\x) = T7^EK^(-)I^>| 2 ( 4 °) 
(-) s 

with the normalization coefficient NTs given in the equation (flop below. 

We can now work out explicitly the probability amplitude for a particle initially in the 
state ip to arrive at x > a in a time t. It will be given by: 



(txs\V r H \ip) = J dE (txs\Er(-))(Er(-)\ip) (41) 

There are two factors in the integrand of the rhs. The first one is obtained from the Eq. 
(|59| ) of section 2: 

(txs\Er(—)) = —p— eX p {-iEt + ixV2mE) (42) 
v 2ir 



The second factor can be computed using the solution ( |55] ) to the Schrodinger equation: 

(Er(-M) = Jdq (Er(-)\q) (<#> (43) 

We said before that (q\ip) practically vanishes for q > 0. Thence, the main contribution 
to the integral comes from q < 0, where (q\Er(— )) = \ frnj 2np T(p)e tpq , with p = \ZlmE. 
In these conditions we obtain 

(Er(-M~JjT*(p)i;(p) (44) 



where ip{p) is the Fourier transform of (q\ip)- Finally, the sought for probability amplitude 
is 

6 r I in ~ 
{txs\V r { ^\i>) ~ J dE J— T *(P) $(p) exp{-iEt + ipx) (45) 

Observe that (txl\Vf_\\ip) = as expected on physical grounds. The amplitude is not yet 
normalized, and has to be divided by the corresponding normalization factor NT_\ . In fact 

N i-) 2 =Ej dt \(txs\r r H m 2 * j o dp \T( P )\ 2 m P )\ 2 (46) 

i.e. N(_ 2 is the probability that the particle be transmitted, as corresponds to the case 
of conditional probabilities we are considering. If the Gaussian chosen as initial state 
is narrow enough in momentum space, and po not too close to pv = V2mV, one can 
approximate it further to the transmission probability 

N{_ } 2 ~ |T(p )| 2 J o d P m P )\ 2 ^\T(p )\ 2 (47) 
12 



as if all the transmission were at the mean momentum p Q . For particles initially (t = 0) 
in the state ip, the average time of arrival at a point x at the other side of the barrier is: 

| r+oo 

W) = TTT^E / dUKtxslV^m 2 (48) 

Any quantity related to the probability distribution PT_Jt\x) can be easily obtained using 
fl45| ) and (|46| ) or fl4"?|) . Observe that the restriction to gaussian wave packets is unnecessary 
for the validity of (^) and ([RJ). It is enough that ip(q) fulfills the conditions given after 

In Fig. 2 we present the most probable toa for a selected range of barriers of various 
heights and sizes. To focus on the effect of the barriers, we keep the same initial wave packet 
and detector position in all the cases. What we obtain can be summarized as retardation 
for E > V, and advancement for E < V. We observe a close agreement between our results 
and the phase time in the cases shown. This is because the stationary phase method 
(that yields r^) is a good mathematical approximation to the true probability distribution 
in these cases in which arg (T(p)) varies slowly. We have also investigated numerically (but 
do not show here) the cases E ~ V where the validity of the approximation is uncertain. 
In fact, the values of obtained in these cases run away from the true maximum of 
PT_s(t\x), which turns out to be very sensitive to the parameters (p Q ,5,V,a) involved. In 
the cases shown, the maxima give excellent approximations to the mean toa t x (i/j), even if 
the coincidence is not exact due to a small amount of skewness present in the distributions. 

In the upper part of Fig. 3, we show the probability distributions in toa P^(t\x) for 
three cases: free propagation, a low V barrier, and transmission by quantum tunneling 
E < V. They look quite similar, apart from the retardation or advancement pointed out in 
Fig. 2. In fact, almost all of the many cases we have investigated have a similar Gaussian 
look. It is only in the cases in which pq5 » 1 does not hold that new features may arise. 
But these are very unrealistic one-particle wave packets, and the results obtained from 
them unsound. We have analyzed in detail the mechanism of advancement in quantum 
tunneling. Usually it is argued that the higher momentum components are transmitted 
preferably through the barrier, and that they propagate faster than the lower momentum 
parts of the packet, being this the cause of advancement. The argument is untenable for 
two reasons: First, advancement is obtained even in the cases of packets whose spread is 
within a range of momenta for which the transmission probabilities |T(p)| 2 are constant 
(i.e. no preferred transmission for higher p's). Second, would the higher momenta be 
transmitted preferably, then the transmitted average momentum should be shifted up- 
wards, and then, the transmitted packet should move faster than the incident one! Of 
course this is physically untenable and false too. 

We now turn to the time of arrival of states reflected from the barrier. First of all, 
we place the detector at a position y < to the left of the barrier. We keep the same 
barrier and the same initial state. Now, it is V l { _ } = I - = J dE U(El(-)) that 
projects on the subspace of the reflected states in which the particle is detected at y in 
some time t. The corresponding probability amplitude can be given borrowing Eq.(f41~1) 
from the transmission case, 

(ty 8 \P l H \iP) = JdE (tys\El(-))(El(-M) (49) 
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There are again two factors on the right hand side. The first one is given by the formalism 
of section 2 as in (£2). Note that, according to fl59p, there will be a change in sign in the 
second term of the exponent, as now r — > I. The second factor , written as before, is 



(El(-m = J dq (El(-)\q) (50) 

In the relevant region q < 0, one has from (|35| ) that (q\El(— )) = ^Jm/2iip (e~ ipg +R(p)e ipq ), 

where p = \j2mE. We now obtain the sum of two contributions (coming from 1 and i?), 
giving 

(El(-M ~ ] j^^(-p) + R*(p)^p)) (51) 

As there is no way that a detector can discriminate between these two contributions, they 
add in the amplitude, and may produce potentially detectable interferences. However, for 
a properly chosen initial wave packet with positive mean momentum such that poS » 1 
or, in other words, with negligible negative momentum tail, the first term in the sum on 
the rhs will practically vanish. Dropping this term, we get the probability amplitude as 
given by: 

(tys\V l H \iP) -^LjdE^ R*(p) i>{p) exp(-iEt - ipy) (52) 

which is a look-alike of the transmission case ([45]). The normalization factor Nl\ is also 
similar: 

N (-) 2 = Y,J dt \(t,y,s\V{^)\ 2 ~ fdp \R(p)\ 2 \i,(p)\ 2 ~ \R(p )\ 2 (53) 

where the last approximation is valid only when the initial state ip{q) is narrow enough 
and, as for transmission, p not too close to py. We can now give the average time of 
arrival at a point y to the left of the barrier of the outgoing particles that were in the 
initial state ip at t — 0: 

/+oo "j r+oo 

dttP^y) = — -Y, / dttKtyslV^m 2 (54) 

with a probability distribution P}Jt\y) that is also normalized to one. In the lower part of 
Fig. 3 we show the probability distributions for reflection by the same potential barriers 
of the upper part. When the barrier acts as a classical mirror (the long dashes curve) 
the distribution is almost the same as that of a free particle that would travel the same 
distance. However, dramatic differences appear for purely quantum reflection (i.e. when 
E > V), whenever the condition p' a = nir, n = 1,2,... (where p' = Jp^ — k% is the 
momentum inside the barrier) holds. In these two bump structure with a dip in 

the middle appears in the distribution. The reason of this behaviour can be found in the 
value of the reflection coefficient 

R(p) = = S ' = ^ exp(z7r — i arctan( — ^— - tan(p'a))) (55) 

^2p 2 p' 2 + k*sm 2 (p'a) 4 PP 
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which vanishes whenever the above condition is satisfied. The proportionality between the 
reflection probability and sin 2 (p'a), that originates the dip, could also be used to explain 
the bumps as due to reflection of the components of momenta greater and lesser than 
p . Were this the case, the bumps should be reflected with different velocities, one above 
and the other below po/m\ Of course, this does not hold. The actual mechanism of the 
bumps is the spread of the incident wave packet in configuration space (not in momentum 
space). The centre of the packet (placed initially at q — go) is reflected with probability 
zero, while the bumps are composed by the parts ahead and behind it (both carrying 
main momentum po). Here, as in the case of tunneling, there is a matching between the 
real and imaginary exponents in the probability amplitude. Thus, while the participating 
momenta are in the small range po ± 8 , due to the Gaussian shape of the wave packet, 
the distances and times involved are macroscopic. Therefore, the phases of the different 
contributions within this range of p are rapidly varying and their sum cancels out. The 
only remaining survivors come from the stationary phases p ~ m(x — go — ar g'(^(p)))/^ 
or p ~ m{—y — g — axg' (R(p))) / 1 . These contributions have still to be weighted by 
the different remaining factors, thus giving rise to the observed probability distributions, 
in which p is a function of t, and depends parametrically on x, (y) and go according to 
the above relations. To illustrate the result of this mechanism, it is enough to point 



out the approximate relation P l ^(t\y) oc sin 2 (yp (^) 2 — ^ 2 °) exp(— 28 2 (p (t) — Po) 2 ) with 
Po(t) = m(—y — go — arg' R(p Q ))/t, as the main responsible of the two bump structure. 

We now recall that the toa of incoming states can be described with the same for- 
malism used above. It is now necessary to project on the subspace of incoming states 
using U(Es(+)). The physical picture here is a detector placed before the barrier to in- 
tercept these states. There are two opposite incoming directions. Mathematically, they 
are selected by the two projectors U(Er(+)) (for detecting right movers at the left of the 
barrier), or U(El(+)) (for left movers at the right). We choose as an example the first of 
these cases. The probability amplitude in times of arrival at y of the incoming states is: 

(tys\V r {+) m = JdE {tys\Er{+)){Er{+m (56) 
Now, instead of ([51]), we have 

(Er(+M ~ l-^(p) + R*( p )^(-p)) (57) 

As occurred for reflection, there are two potentially interfering terms, due to the lack 
of discrimination between the direct incoming component and the receding part of the 
reflected component. However, only the first term will give a sizeable contribution when 
the negative momentum tails are negligible. It is worth noting that the direct component 
gives the same time of arrival distribution that in the absence of barrier, as it should be 
on physical grounds. In the stationary phase approximation, the peak of the distribution 
is on the classical toa: t = m(y — qo)/po, as expected. This also gives the meaning of 



the negative toa's [Q: they are obtained when the detector is placed to the left of the 
incoming right mover, so that the particle would have arrived at it before t = 0. Finally, 
we note that the probability amplitude of incoming states at a point x to the right of the 
barrier would be proportional to the receding part of the transmitted state. 
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We now restore the (±) labeling that should be explicitly included in the toa formu- 
lation of Section 2, and that we ignored from Eq. flI%D on. By doing it, instead of the 
equations (^7|,^9|) at the end of that section we would arrive to 



, r+co 

t x (±) = J2 dtt\txs(±))(txs(±)\, (58) 

with the position dependent eigenstates given by 

\txs(±)) = dE exptiEt - isxy/2mE)\Es(±)). (59) 

\/2tt Jo 



. 2tt Jo 

We thus see that many of the physical considerations made in this section in order to 
obtain the probability distributions in toa are automatically taken into account by the 
formalism. Instead of Eq. ( fiO|) we would get directly 



1 

'(-) 



P(^(t\x) = \(txr(-m\ 2 (60) 



and similar abridged (but fully equivalent) expressions in analogous cases. We recall from 
a formal point of view, that there are two possible constructions of £(±) starting from the 
free toa, as there are two constructions of the stationary scattering states \Es(±)) starting 
from the free states \EsO). Even if each of these, £(+) or t(— ), is associated to a kind of 
asymptotic state, in or out respectively, nowhere in the construction have we imposed or 
used any asymptotic value for x, y, t or for any other quantity in the problem. It is the 
finite extension of the barrier the only responsible for the simple asymptotic looking of our 
expressions. We used this to discuss the physical questions with a minimum of burden. 



5 Conclusions 

We have worked out a formalism for obtaining the time of arrival at a space point of par- 
ticles that move through interacting environments. Our construction follows a circuitous 
path; we desist from first computing the classical toa of the problem, and then quantize 
it, a procedure that leads to a dead end. Instead, we start from the quantum toa of the 
free moving particle, and then transform it canonically to the interacting case. This is 
achieved by the use of the appropriate Moller operators, that we employ in a nonstandard, 
but orthodox, setting. Usually, they are introduced to connect the values (at t = 0) of 
those free and scattering states that will coincide asymptotically (at t — > ±oo, depending 
on the case). Here, we first derive them from the connection between the free and inter- 
acting Hamiltonians, and then use them to connect the free and interacting fcoa's, which 
solves our problem. For simplicity, we have only addressed explicitly the cases in which 
the transformations are unitary. The well known results from formal scattering theory 
pave the way to the general cases that require isometric transformations. 

We have also performed a quite exhaustive (and a bit heavy) analysis of the relation 
between classical trajectories and scattering states. We found necessary to delve into 
the parallelism between the partitions of classical phase space by means of in or out 
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trajectories on the one side, and the spans of the Hilbert space by means of the (+) or 
(— ) states on the other. The right appreciation of the different issues involved in that 
correspondence was necessary for a successful use of the formalism previously developed. 
At the end, we reached the conclusion that this formalism did incorporate from the outset 
these subtleties. In fact, for the asymptotically vanishing potentials considered here we 
obtain two different operators t(+) and t(— ). They measure the toa at a space point of 
the incoming and outgoing states respectively. 

In the course of our numerical analysis we have detected that the phase time t$ gives 
a good approximation to the most probable time of arrival. It provides a first estimate 
of the time spent in the transmission or reflection, after subtracting the time of free 
flight. We have found advancement of the transmitted wave packet in the case of pure 



quantum tunneling. This is a long known phenomenon first predicted by Hartman ||19| , 
and experimentally evinced by the two photon experiments at Berkeley [|17|, |18| . We have 
found an unexpected phenomenon for purely quantum reflection: the two bump structure 
that appears when p' Q a = rnr. We have shown it only in one case in Fig. 3, which by 
no means is the most crisp case. Neat double humps appear in all the expected cases, 
the lesser the potential barrier the finer the feature. With the insight provided by our 
formalism, we have also proposed an explanation to this structure. We think this feature, 
even if less spectacular than the superluminal tunneling of photons, deserves experimental 
confirmation. A minor modification of the two photon experiment could serve for this 
purpose. It would be enough to place a quantum mirror in the path of one of the entangled 
photons, and check for the presence (or absence) of a two dip structure in the number of 
coincidence counts. 
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p 




Figure 1: Alternative partitions of the phase space plane (q,p) in terms of the trajecto- 
ries of particles with Hamiltonian H(q,p) = p 2 /2m + V/ cosh 2 (g/<f). The trajectories of 
incoming right movers r(+) are above the separatrix p = pv(q), and those of incoming 
left movers below it, both in the upper figure. The trajectories of outgoing right 
movers r(— ), and left movers /(— ) are respectively above and below p = —pv(q) in the 
lower figure. The notation is borrowed from that of the stationary scattering states 
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Figure 2: Most probable time of arrival and phase times for transmission through barriers 
of different heights and widths. Times are represented on the vertical axis against the 
barrier width ranging from to 10. The initial wave packet (c.f. Eq. fl3~9"|)) has go — 
—50, po = 2, 5 = 10 and m = 1. The detector is at x = 50. These parameters are fixed, 
so that the classical toa for free propagation is 50 in all cases, (with our packet, we get the 
most probable toa at to = 49.9377 quantum mechanically). We show our results for the 
most probable toa with triangles, squares, pentagons, and crisscrosses, corresponding to 
barrier heights of V — 0.5, 1.125, 3.125, and 4.5 respectively. The curves represent "to+ 
Wigner time delay" for these same heights, (being the solid one for V = 4.5, etc). The plots 
neatly indicate retardation when transmission is classically possible, and advancement for 
pure quantum tunneling. 
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Figure 3: Probability distributions in times of arrival for transmission (above) and re- 
flection (below), corresponding to equations ( |40"1) and (|54]) . The initial wave packet is in 
all cases the same of Fig. ([$]). The total length of the particle's path is always 100 units, 
so that the classical free toa is 50. The width of the barrier is a = 10. The solid line 
represents free transmission, the short dashes the cases with V = 1.125, and the long ones 
those with V = 4.5. The two bump structure is a general result in the cases where there 

is no classical reflection as explained in the text. 
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